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The rotated multiplier method is performed in the case of the boundary stabilization of the 
wave equation by means of a (linear or non-linear) Neumann feedback. This method leads to new 
geometrical cases concerning the "active" part of the boundary where the feedback is applied. Due to 
mixed boundary conditions, these cases generate singularities. Under a simple geometrical condition 
concerning the orientation of the boundary, we obtain a stabilization result in both cases. 

Introduction 



' In this paper we are concerned with the stabilization of the wave equation in a multi-dimensional body 

by using a feedback law given by some part of the boundary of the spacial domain and some function 
defined on this part. The problem can be written as follows 



in Slxi;, 



u = Q ondSl D xE* + , 

d v u = F on <9f2jv x K!J_ , 

it(0) = uq in Q , 

u'(0) = u\ in Q , 



where we denote by u', u" , Am and d v u the first time-derivative of u, the second time-derivative of the 
scalar function u, the standard Laplacian of u and the normal outward derivative of u on dfl, respectively; 
t^. \ {9Q.D, dftjj) is a partition of <90 and F is the feedback function which may depend on the state(it, u'), 

■ the position x and time t. 

Our purpose here is to choose the feedback law, that is to say the feedback function F and the "active" 
part of the boundary, BQn, so that for every initial data, the energy function 



E{u,t) = \ f {\u'(t)\ 2 + \Wu(t)\ 2 )dx : 



is decreasing with respect to time t, and vanishes as t — > oo. 
Formally, we can write the time derivative of E as follows 



E'(u,t) = / Fu da, 



and E is non- increasing if Fu' < on dfl^- 

Thanks to the multiplier method introduced by L.F. Ho [12] in the framework of Hilbert Uniqueness 
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Stabilization and rotated multiplier method 



Method [llj . it can be shown that the energy function is uniformly decreasing as time t tends to oo, 
under suitable assumptions. 

In this method, one chooses some vector field x i— > m(x) = x — xq and 

dil N = {x E dn / m(x).v{x) > } , F = -{m.v)v! , 
where v is the normal unit vector pointing outward of f2. 

This method has been performed by many authors, especially by Komornik and Zuazua |10| . Following [7] , 
Bey, Loheac and Moussaoui obtained in U similar results under weaker geometrical conditions concerning 
the interface T = dflN H d£ln- We here extend the above result for rotated multipliers introduced by 
Osses PUCE]. 



1 Notations and main results 

Let f2 be a bounded open connected set of M. n (n > 2) such that 

dQ is C 2 in the sense of Necas [H] . 



(1) 



Let xq be a fixed point in R™. We denote by / the n x n identity matrix, by A a real nxn skew-symmetric 
matrix and by d a positive real number such that d 2 + \\A\\ 2 = 1. We now define the following vector 
function: 

Viet", m(x) = (dl + A)(x-x ). 
We consider a partition (cTJjv, dtlo) of dfl such that 



(2) 



r = 8Qd n 9fijv is a C -manifold of dimension n — 2, 

m.v — on F, 
dfl flu is a C 3 -manifold of dimension n — 1, 
H n - l {dSl D ) > 0, 

where cu is a suitable neighborhood of V and H n ~ 1 denotes the usual n — 1-dimensional Hausdorff measure. 
Let g : K — > R be a measurable function such that 



g is non-decreasing, and 3K > 0; \g(s)\ < a.e. 
Let us now consider the following wave problem 



(3) 



u" - Au = 
it = 

(S) I d v u = -{m.v)g(u') 
u(0) = u Q 

^ U (0) = Ml 



n x 



on x R* + , 
on OfiAT x m 
in Q, 
in fi. 



for some initial data 



(uo,Ui) 6Hj,(!!) xL 2 (!l) := H 
with Hj,(n) = {ue H x (f2) ; w = on 30^}. 

This problem is well-posed in this space. Indeed, following Komornik 9J, we define the non-linear operator 
A on H by 

A(u,v) = (—v,—Au), 

D(A) = {(u,v) £ H^(O) x H^(O) ; Au e L 2 (fi) and = -(m.v)ff(w) on 5^}, 



so that (S) can be written in the form 



(u, v)' + A(u, v) = 0, 
(u,v)(0) = (u ,ui) . 



It is classical that A is a maximal- monotone operator on H and that D(A) is dense in H for the usual 
norm. Following Brczis [Q, we can deduce that for any initial data (uq,vq) G D{A) there is a unique 
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strong solution (u, v) such that u G W 1 '°°(R; Hjj(fi)) and Au G L°°(R + ; L 2 (S1)). Moreover, for two initial 
data, the corresponding solutions satisfy: 

vt > o, iK^t),* 1 ®) - («*(*), « 2 (f))||* < ll(4^o) - ■ 

Using the density of D (A), one can extend the map: 

(u o ,v ) ' — > (u(t),v(t)) 
D(A) — > H 

to a strongly continuous semi-group of contractions (S(t)) t >o and define for (uo, vq) G H the weak solution 
(u(t),vf(t)) = S(t)(u ,ui) with the regularity u G C(M + ;Hd(^)) n C 1 (M+; L 2 (^)). We hence define the 
energy function of solutions by: 

E(u,0) = ^ [ (|ui| 2 + \Vu \ 2 )dx, E{u,t) = \ [ {\u'(t)\ 2 + \Vu(t)\ 2 )dx if t > 0. 

^ JQ. & JO. 

In order to get stabilization results, we need further assumptions concerning the feedback function g 

3p>l,3/c>0, \g{s)\ > fcmin{|s|,|s| p }, a.e. . (4) 
Concerning the boundary we assume 

dn N c {x G dCl, m(x).u(x) > 0}, dn D c {x G 917, m(x).j/(x) < 0}, (5) 
and the additional geometric assumption 

m.T < on r , (6) 

where t(x) is the normal unit vector pointing outward of 8Qn at a point x € T when considering 8Qn 
as a sub-manifold of dfl. 

Remark 1 It is worth observing that it is not necessary to assume that 

H n ~\{x G dn N ; m{x).u{x) > 0}) > 0. 

to get stabilization. In fact, our choices of m imply such properties ( see examples in Section 4 ) whether 
the energy tends to zero. 

Following |12j . Rellich type relation [17] is useful in control of the wave problem. Then, Komornik and 
Zuazua [lOj have shown how these relations can also help us to stabilize the wave problem, but failed to 
generalize it in higher dimension than 3. The key-problem was to show the existence of a decomposition 
in regular and singular part ( [BJ [H] ) which can apply to stabilization problems or control problems in any 
dimension. The first results towards this direction are due to Moussaoui [13] , and Bey-Loheac-Moussaoui 
[1] who also have established a Rellich type relation in any dimension. In this new case of Neumann 
feedback introduced by Osses [T51 [TB], our goal is then to generalize those Rellich's relations to get 
stabilization results about (S). This will lead us to get a stabilization result under (O, ©. 
As well as in [3] , we shall prove here two results of uniform boundary stabilization. 

Exponential boundary stabilization 

We here consider the case when p = 1 in This is satisfied when g is linear: 

3a > : Vs G E , g(s) = as . 
In these cases, the energy function is exponentially decreasing. 

Theorem 1 Assume that geometrical conditions (J2]), © hold and that the feedback function g satisfies 
© and (g|) with p = l. 

Then under the further geometrical assumption © , there exist C > and T > ( independent of d) such 
that for all initial data in H, the energy of the solution u of satisfies 

Vt>^, E(u,t) < E(u,0) exp(l- ^tj . 
The above constants C and T depend only on the geometry. 
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Rational boundary stabilization 

We here consider the general case and we get rational boundary stabilization. 

Theorem 2 Assume that geometrical conditions ((2|), © hold and that the feedback function g satisfies 
© and dU) withp > 1. 

Then under the further geometrical assumption ^ , there exist C > and T > ( independent of d) such 
that for all initial data in H, the energy of the solution u of satisfies 

Vt>^, E(u,t) < G . 
where C depends on the initial energy E(u,0). 

Remark 2 Taking advantage of the works of Banasiak- Roach J^j who generalized Grisvard's results 
in the piecewise regular case, we will see that Theorems]^ and\^ remain true in the bi- dimensional case 
when assumption ([I]) is replaced by following one 



dfl is a curvilinear polygon of class C 2 , 
each component of dQ \ T is a C 2 -manifold of dimension 1 



(7) 



and when condition © is replaced by 

ViGT, < vox < 7r and, if w x = 7T, m(x).T(x) < . (8) 
where vo x is the angle at the boundary in the point x. 

These two results are obtained by estimating some integral of the energy function as well as in [S] . This 
specific estimates are obtained thanks to an adapted Rellich relation. 

Hence, this paper is composed of two sections. In the first one we build convenient Rellich relations and 
in the second one we use it to prove Theorems [T] and [H 

2 Rellich's relation 

2.1 A regular case 

We can easily build a Rellich relation corresponding to the above vector field m when considered functions 
are smooth enough. 

Proposition 3 Assume that f2 is an open set ofW 1 with boundary of class C 2 in the sense of Necas. If 
u belongs to H 2 (f2) then 

2 / Au (m.Vu) dx = d(n - 2) / \Vu\ 2 dx+ I {2d v u (m.Vu) - (m.v) |Vu| 2 ) da . 
Jn Jn Jan 

Proof. Using Green-Riemann identity we get 

2 / Au{m.Vu)dx = I 2d v u{m. Vu) da - 2 / Vu.V(m.Vu) dx. 

Jn e JdQ. e Jn e 

So, observing that 

Vu.V(ra.Vu)) = (D 2 u)(m, Vu) + {dVu + AVu).Vu = im.V(|Vu| 2 ) + d\Vu\ 2 + AVu.Vu 

{D 2 u is the Hessian matrix of u) and since A is skew-symmetric, we get 

2 / Au(m.Vu)dx= / 2d v u{m.Vu) da - 2d / \Vu\ 2 dx - / m.V(| Vu| 2 ) dx. 
Jn e Jon,, Jn e Jn e 

With another use of Green-Riemann formula, we obtain the required result for div(m) = nd. ■ 

We will now try to extend this result to the case of an element u belonging less regular when is smooth 

enough. 
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2.2 Bi-dimensional case 

We begin by the plane case. It is the simplest case from the point of view of singularity theory and its 
understanding dates from Shamir [IS]. 

Theorem 4 Assume n = 2. Under geometrical conditions ^) and ([7]), let u G H 1 (J7) such that 

AueL 2 (f>), u /SilD eH 3 / 2 (a!! D ), d v u /dQN eR^ 2 (dn N ). (9) 

Then 29„u(m.Vu) — (m.^)|Vu| 2 6 L 1 (057) and there exist some coefficients (c x ) x ^r such that 

2 Au(m.Vu) dx = / (29„u(m.Vu) - (m.t/)|Vu| 2 ) der + - V] c 2 (to.t)(:z;). 
Jan 4 

Proof. We first begin by some general considerations which will be used in the general case too. It is a 
classical result that u G H 2 (w) for every open domain to such that u> (s Cl\T. For sake of completeness, 
let us recall the proof. 

A trace result shows that there exists ur G H 2 (uj) such that ur = u on cT^d and (9„u.r = <9„u on OSIn- 
Hence, setting / = Amjj — Am G L 2 (r2), u$ = u — ur satisfies 

—Aus = f in SI, 

u s = ondQ D , (10) 

<9„m,s = on dflN- 

Now, if w <e fi\T U c?Slz3 and £ is a cut-off function such that £ = 1 on to and supp(£) C SI, then, for a 
suitable g G L 2 (S1), it w = us£ satisfies the Dirichlet problem 

Au u = g on f2, 
u w = on <9S1, 

and using classical method of difference quotients ([H]), one can now conclude that u w G H 2 (S1), hence 
«5 G H 2 (w). Else, if w <g Sl\r U 90jv, and £ is a cut-off function such that £ = 1 on to and supp(£) C fi, 
then, for a suitable g G L 2 (S1), it w = ug£ satisfies the Neumann problem 

—Auu + u w = g on SI, 
d„u w = on <9S1, 

and, using similar argument, one gets us G H 2 (w). 

Let Sl e = {x G SI; d(x, T) > e}. By compactness of Sl £ , we get u G H 2 (S1 £ ). An application of Proposition 
131 to our particular situation gives us the following relation 

2 / Au(m.Vii) dx = / (2d v u(m. Vu) — (m.i/)| Vu| 2 ) da 

and we will try to let e — > 0. Using derivative with respect to and t, we get 



Au(m.Vu) dx = / (m.v) {{d u u) 2 — (d T u) 2 ) da + 2 {m.T){d v u)(d T u) da. 
Jan e Jd(l e 

First, since Am G L 2 (S1) and u G H^Sl), Lebesgue's dominated convergence theorem immediately gives 

lim / Au(m.Vit) dx = I Au(m.Vu) dx 
Jn E Jn 

Now, we work on boundary terms. We shall introduce the following partition of 9Sl e : let <9Sl e = <9Sl e n9Sl, 
<9S1* = dfl e (~1 SI . and a decomposition result. 

We use a decomposition result due to Banasiak and Roach [2]: every variational solution of (fT0|) can be 
split as a sum of singular functions. There exists some coefficients (c x ) xe r and ur G H 2 (f2) such that 

u = u R + ^c x U§ =:u R + u s (11) 
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where Ug are singular functions which, in some neighborhood of x G T, are defined in local polar 
coordinates by: 

U§(r,0) = p(r)r^ sin(- 6). 

LVJ X 

with p some cut-off function(sec Fig. 1). 

Using the density of C 1 (fi) in H 2 (17), we will be able to assume that u R G C 1 (r2). Let us look at boundary 
terms on dVl e first. We claim first that for some constant C > , 

\m.u\ < Cd(.,r). 

In fact, if x G f2 and x\ G T which satisfies \x— x\ \ = d(x,T), one gets: 

m.i>(x) = m(x).(i>(x) — v{x\)) + (m(x) — m(x\)).v(xi) (observe that m.v(x\) = 0). 
Hence, using the fact that v is a piecewise C 1 function (see Fig. 2), we get 

\m.v(x)\ < (HmlloolKHoo + l)d(a:,r). 
Now, working in local coordinates, one gets 

d(x,T)\Wu\ 2 G L°°(0fi). 

Hence Lebesgue theorem implies 

lim [{m.v)((d u u) 2 -(d T u) 2 )da = [ (m.v) {{d u uf - (d T u) 2 ) da. 



on e Jan 



On the other hand, assumptions ([9]) give d v u G H 1//2 (9J7n), d T u G H 1 / 2 (<90n) and 
d v u G H- 1 /2(^ D ) i q tU g H 1 / 2 (5r2 D ). Hence we get 



{m.T)(d u u){d T u)da — > / {m.T){d v u){d T u)da. 

on e Jan 

Now, we have to consider the boundary term on <9f2*, 7 e (Vu). This quadratic form in Vit can be 
decomposed using pip in 

I e {Vu R ) + 2J e {Vu R , Wu s ) + I e (Vu 8 ), 

where J e is the corresponding bilinear form. 
Concerning I £ (\7u R ), regularity of m gives the estimate 

\Ie(Vu R )\<C [ \Vu R \ 2 da 

J dVt* 

This term is then an 0(e) since Vu R is bounded on fl. 

For the term I E (\7us), we first note that, adjusting the cut-off functions, the supports of u x s and u v s are 
disjoint. Hence, using decomposition (fTTj) . we can write 

I e {Vu S ) = Vc 2 f (2dvU x s (m.\7u x s ) — (m.v)\Vu x s \ 2 )da. 

If vo x < 7r, one has 

2d v u s {m.Vu s ) - {m.v)\Vu x s \ 2 = 0(e^ 2 ) on C e {x). 

Hence, after integrating on C e (x), we get lim e ^ If (e) = 0. 
If w x = 7r, we will need the following identity 

2d v u x s (m. V115) - (m.j/)|Vus| 2 = -^-to.t(:e) on C e (x). 

If we observe that C e (x) behaves as a half-circle when e — > 0, an integration gives 

lim / (2(i/. Vit|) (m. Vu|) - (TO.i/)|Vus| 2 )dcr = jm(x).r(x). 



C e {x) 
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Figure 1: Shape of the boundary near an angular point x. 

T (y) 

dn D 




Figure 2: Unit vectors v{x), t(x), v{y) and r{y) when dVL is regular at x. 

Finally, the bilinear term J £ (\7ur, Vits) can be written entirely 

/ d„un(m.X7us) da + / d v us{m.Vun) da — I (m.i')VuR.'Vus da. 
Jan* JdO-i Jan* 

Using the regularity of m and Cauchy-Schwarz inequality, we get an estimate of the form 



\J e (Vu R ,Vu s )\ <C( / \Vu R \ 2 da 
'on* 



1/2 



( [ \Vu S \ 



da 



1/2 



We have seen that the first term in this inequality vanishes when e — > 0. For the other term, we now 
observe that, if e is small enough 

on; = y c e (x) 

where C E {x) is an arc of circle of radius e centered in x. Then, we may write 

/ \Vu s \ 2 da<2 V c 2 y [ \VU v s \ 2 da. 

A similar computation shows that, for x € T, / \VUg\ 2 da = O(l). Moreover, if x ^ y, f/| is bounded 

JC E (x) 

near x, we get / \VUg\ 2 da = O(e). This completes the proof. ■ 



Remark 3 The assumption 7i l (dVlr>) > is not necessary in the above proof. We will now see why we 
need this assumption on the Dirichlet part in higher dimension. 
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2.3 General case 

We now state the result in general dimension. 

Theorem 5 Assume n > 3. Under geometrical conditions ([T]) and let u G H 1 (J7) such that 

AueL 2 (fi), u /dnD eR 3/2 (dn B ), d v u m „ e R^ 2 (dn N ). (12) 
Then, 2d„u(m.Vu) - (m.^)|Vu| 2 G L 1 (9Q) and i/iere eiiste C G H 1 / 2 (r) sucfc tftai: 



2 / Au(m.\7u)dx = d(n-2) / Vu| 2 efe + / {2d v u(m.\Ju) - (m.z/)| Vw| 2 ) dcr + / (to.t)|C| aV 
./n Jn Jan Jr 

Proof. We will essentially follow [3]. As in the plane case, we let fl e = {x G ft; d(x,T) > e}. For any 
given e > 0, we have the identity of Proposition [3] 

2/ Au(m.Vu) dx = d(n - 2) \Vu\' 2 dx + (2d v u{m.Vu) - (m.u)\ Vu\ 2 ) da, 
Jn E Jn E Jan e 

and we will again analyze the behavior of each term as e — > 0. 

First, since Am G L 2 (ft) and u G H 1 (ri), Lebesgue's dominated convergence immediately gives 

lim / Au(m.Vu) ds = / Au(m.Vit) dx, lim / |Vw| 2 da; = / |Vu| 2 da;. 
£ ^°Jn E Jn £ ^°Jn e Jn 

Next, we will work on boundary terms. 




Figure 3: Picture of dft* and dft e 

We define again <9fL = dfl e (~l dfl and <9f2* = dfl e D ft (see Fig. 3). As well as in the plane case, there 
exists some constant C > such that \m.v\ < Cd(., V). Thus, using the fact that 

d(,r)|Vu| 2 G 1^(00) 

(see [I], Proposition 3), we can use again Lebesgue's theorem to conclude that, as e — > 0, 



(m.v)\Vu\ da — > / (m.v)\S7u\ da. 
an c Jan 

For the second integral, defining Vgsi as the gradient of derivatives along dfl, we write that 

d v u (m.Vu) = {m.v) |<9„w| 2 + d v u (m.Vanit). 

The first term is integrable. The second one is, on 8ft n, the product of a H 1 / 2 term by a H" 1 / 2 one and, 
on dfti), the product of a a H -1 / 2 term by a H 1 / 2 one. Hence, Lebesgue's theorem gives again, as e — > 0, 



d v u(m.Vu)da — » / d„u(m.Vu)da. 
an e Jan 
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Now, wc consider the terms on cWl*. We assume that e ^ £o and we define: ui eo :— f2\f2 £o . Now, as well 
as in the plane case, we can write 

u = u R + u s (13) 

where us is the variational solution of some homogeneous mixed boundary problem and ur belongs 
to H 2 (i*j £o ). Working by approximation if necessary, we can suppose that ur G C 1 (w £o ). Consid- 
ering the same quadratic form as in the bi-dimcnsional case, this leads to the following splitting of 

{2d v u(m.Vu) — (m.v)\ Vu\ 2 )do~: 

Z e (Vujt) + I e (Vu s ) +2J £ (Vu R ,Vu s ). 

Sice Vui? 6 L°°(a; eo ) and TL n ~ l {dW e ) — > as e — ► 0, it is obvious that the first term I £ (\7ur) tends to 
zero. 

The bilinear term J e (Vwfl, Vus) is, as above, 

/ d u un{m.Vus)da + / d^us(m.Vuji)da — / (m.f)Vuji.Vusda. 
Jon* Joa* JdO-i 

Using the regularity of m and Cauchy-Schwarz inequality, we get an estimate of the form 

|J e (V«K,Vus)|<c(/ \Vu R \ 2 da) (/ \Vu s \ 2 do) (14) 

As above, it is clear that the first term tends to zero as e — > 0. In order to analyze I e (Vus) we will need 
further results. 

To begin with, we shall use some definitions. Every x G dfl* belongs to a unique plane x* + (t*,v*) 
(setting: r* = t(x*), v* — i/(x*)) and more precisely to an arc-circle C e (x*) in this plane of center x* G T 
and of radius e (the figure is similar to Fig. 2 in the plane x* + (r*, v*) ). Wc define 

D e (x*) :=Lu e n(x* +(t*,p*)). 

For any x G D eo (x*), we separate the derivatives of u along the sub-manifold x — x*+T with the co-normal 
derivatives: 

Vu(x) = V r u(x) + V 2 u(x), V r u(x)eT x ,T, V 2 u(x) G (r* , i/*) . (15) 

Using methods of difference quotients (see for instance [3], Theorem 4), one gets Vrit G H 1 (w £o ) i.e. 
Vr«s G H 1 (w Eo ). We will also need the following result concerning the behavior of boundary integrals 

Lemma 6 Let Eq > 0. Assume that u is such that u = on duj £o D dflo, 

Vi'er, ii(i , ,.)GH 1 (D £0 (x')) 

and 

(x*^||«(,OIIhi(d. cx«))) eL 2 (r). 

TTien t/iere exists C > depending only on fi sweft i/iaf, /or any e sufficiently small, 

\\ u { x \-)\\h(C e (^)) d l ^ lk(a ; *>-)llHi(D e (x*)) d 7(^*)- 

Proof of Lemma HH We begin by changing coordinates as well as in [4]. For every Xq G T, there exists 
po > 0, a C 2 diffeomorphism from an open neighborhood W of Xq to B(po) ■= B n _2(po) x -^(po) (see 
Fig. 5) such that 

0(4) = o, 
e(Wnfi) = {2/gB(po);2/»>o}, 

e(^n9firj) = {y £ B(po); y n -i > 0, j/„ = 0}, 
e(VFnafijv) = {y G B(po); y n -i < 0, y n = 0}, 

9(wnr) = {ye%);!/„_i = 0,2/„ = 0} := 7 (p ). 
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Figure 4: The set W. 

Reducing e if necessary, we may assume that -D £o (:z;q) C W. We then get, writing for x £ W, Q(x) = 
(y' , y) and v := u o © _1 : 

( f u 2 di] di(x*) = I ( [ v 2 d£(yj) dy' . 

Setting 

BUp) ■= iv = (Vn-uVn) G B 2 {p);y n > 0}, C+(p) := {y = (y n -i,y n ) G dB 2 (p);y n > 0}, 

we first observe that we can choose /o^* such that {?/} x B^ip) C 0(D E (a;*)). Hence if we denote by 7r 2 
the projection on {0^-2 } x R 2 , the change of variables 



7r 2 (6(C £ (x*))) 



c 2 + (p) 



gives the estimate 

f v(y' ,y) 2 d£(y) < C f v(y' , zf di{z) (16) 

Je(c e (x*)) Jc+{p) 

for a constant C depending only on Xq. 

We will now estimate this latter integral in terms of ||V21'||^2(/ 3/ »i. x b+( j0 ))- If v p{y) '■= v(y',y), one gets 



Vv p 6 L 2 (B£(1)) and 



l|V^|| L2(B +(l)) = l|V2«|| L2({y , }xB +( p )), IMl=(C+(1)) = P *W V hmy>}xCt(p))- 

Observing that v p = on i? 2 f+ (l) := {(y n -i,yn) G B^"(l);?/ n > 0}, trace theorem and Poincare's 
inequality give, for some universal constant C > 0, the estimate 



c 2 + (p) 



U 2 ( y ',y)d£(y)<Cp||V 2 ^|| 2 2({y;}xB2+(p)) . 



Hence, thanks to (fll)]) . one gets 



e(c e (z*)) 



w 2 ( y ',yW)<C P .-||V 2t ;|| 2 2({ , y;}xB2+( ^ )r 



Observing that p^* is uniformly 0(e) on W (~1 T and the diffeomorphism 0(x*, .) (see Fig. 6), we can 
conclude that, for some constant C x * depending only on Xq 



'e(C e (x*)) 



v 2 (y,y)de(y) < C x *s\\u(x* , ■) ll H i(e-i({ y <}xB+(p))) ^ C x3 £||u(a;*, .)|||i( De (x*))- 
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Q(x*,.) 




v{x*) 



Figure 5: The C 2 -diffeomorphism <d(x*, .) in the plane x* + (r* ,v* 



Hence, after an integration on W (~l T 

( [ u 2 dt)d 7 (x*) < C x *e ( \\u(x*, OIIhicd.Cx.))^^*)- 



We finally complete the proof by using a partition of unity on the open sets (W x * )x*er- 
End of proof of Lemma\(j[ M 

Let us come back to our problem. Using (|15|) for us, Pythagore's theorem gives 



\Vus\ 2 da = / \Vrus\ 2 da 
an* Jan* 



\\7 2 us\ 2 do-. 



(Kit 



Applying Lcmma[6]to Vrus, we get that the first term vanishes as e — > 0. As well as in the bi-dimensional 
case, we will see that the second term above is bounded, using more information on us- 
Thanks to [I] (Theorem 4) and Borel-Lebesgue theorem, we may write 



u s (x) = r]{x*)U s (x - x*) := v(x*)Us (x) on uj e 



(17) 



with Us locally diffcomorphic to Shamir's function, and rj £ i? 1//2 (r). We then get, thanks to Fubini 
theorem 



\V 2 u s \ z da= / t?(x*) 2 / V 2 t/f \ z dl )dj(x*), 
an* Jr K Jc e (x*) ' 

and, as well as in the bi-dimensional case, we show that this term is bounded by O(l) |Mlx,2(rv We have 
now proven that the second term in (|14|) is bounded, that is 

J e (Vu_r) — > as e — > 0. 

To treat the last term J e (Vus), we will use the same tools as above. The splitting (fT5| for us gives us 

4(Vu s ) = / £ (V 2 U5) + J e (V r «s) + 2J e (V 2 M S , Vrus)- 

As above, the term I £ (Vrus) is estimated by / |Vrus| 2 G?cr; it then vanishes for e — > 0. The bilinear 

Jon* 

term is estimated by 

([ \V 2 u s \ 2 da) 1/2 ( f \Vru s \ 2 do-) 

Van* ' van* ' 

it then tends to zero since the first term is bounded and the second one vanishes for e — > 0. 
For the last term I e (V 2 us), we use (fT7|) and Fubini's theorem to write it 

V(x*) 2 ([ , 2(v.\7 2 Uf)(m.\7 2 uf )-( m .v)\\7 2 uf\ 2 )de)dj(x*). 

r Vc £ ((i-) 7 

We first work in the plane a;* + (t* , — v*) and, as above, we get 



1/2 



lim 

e— 



(2(i/.V 2 l/|*)(m.V2?7I*) - (m.i/)|V 2 J7I*| 2 )d« = '-m{x*).T{x*) 
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Moreover, for any e > 0, this integral term on C e (x*) is dominated by ?||m||oo G ^(T); so dominated 
convergence theorem applies and finally 

lim I E (V2Us) — — / rfm.rd'y. 
The proof is now complete with £ = rj. ■ 

We will now apply Rcllich's relation to the stabilization of solutions of (S). 



3 Proof of linear and non-linear stabilization 

We begin by writing the following consequence of Section 2. 

Corollary 7 Assume that 1 1— » (u(t),u'(t)) is a strong solution of (S) and that the geometrical additional 
assumption ^ if n > 3 or ([5]) if n — 2 holds, then, for every time t, u(t) satisfies 

2 / Au(ra.Vit) dx < d(n - 2) / \S7u\ 2 dx + / {2d u u{m.Vu) - (m.u)\ Vm| 2 ) da. 
Jn Jn Jan 

Proof. Indeed, under theses hypotheses, for each time t, (u(t),u'(t)) G D(A) so that u{t) satisfies (O or 
(fT"2ll . The corollary is then an application of Theorem |4] or [5] ■ 



p — 1 

We will be able to prove Theorems [T] and [3] showing that, for a = — - — , one can apply the following 
result ([!]): 

Proposition 8 Let E : R + — > M + a non-increasing function such that there exists a > and C > 
which fulfills 

/oo 
£ Q+1 (s)ds < OE(i). 

Then, setting T = CE a (0), one gets 

ifa = 0, Vi > T, E(t) < E(0) exp ( 1 - - 



i/a>0, Vi>T, E{t)<E{0) 



We come back to our proof now. 

Proof. Following [5] and [5], we will prove the estimates for (uo,ui) £ D{A) which, using den- 
sity of the domain, will be sufficient to get the result for all solutions. 
Setting Mu = 2mSJu + d(n — l)u, we prove the following result. 

Lemma 9 For any < S < T < oo, one gets 

2d f E^dt < - E^ [ u'Mudx + ^—^ [ E^ E' ( [ u'Mudx) dt 
Js l Jn J s 2 j g \J n ) 



+ / E— 

Is 



{m.v)((u') 2 - \Vu\ 2 - g(u')Mu) da) dt. 

J 



Proof of Lemma[9j Using the fact that u is solution of (S) and noting that u" Mu = (u' Mu)' — u'Mu', 
an integration by parts gives: 

E 2 ^ (u" - Au)Mudxdt 



s 



E P ^ I u'Mudx 
o 



^ / E^E'f f u'Mudx) dt- [ E^ [ (u Mu + AuMu) dx dt. 
s 2 J s \Jq ) J s Jn 
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Corollary [7] now gives 

/ AuMudx < d(n-l) I Auudx + d(n-2) / \Vu\ 2 dx + / (2d„u(m.Vu) - ( m.v)\ Vw| 2 ) da. 
Jn. Jn Jn Jan 

hence, Grcen-Riemann formula leads to 

/ AuMudx < -d / \Wu\ 2 dx + / (d v uMu-{ m.^)|Vu| ) da. 
Jn Jn Jan 

Using boundary conditions and the fact that Vu = d v uv on dflo, we then get 



AuMudx < -d \Vu\ 2 dx- (m.v)(g(u')Mu + \Vu\ 2 ) da. 
n Jn Jan N 

On the other hand, using div(m) = nd, an other use of Green formula gives us 

u'Mu' dx = -d / (u') 2 dx+ / {m.v)\u'\ 2 da. 
n Jn Jan N 

End of proof of Lemma\j^ ■ 

Coming back to our problem, Young inequality gives 



u'Mudx\ < CE{t). 

n 

Lemma H] shows that 

2d I E^dt < C{E S ^{T)+E £ ^-{S)) + C l E^E'dt 
Js Js 

,'\2 



■ P-i 



+ / E— / {m.v)((u') 2 -\S!u\ 2 - g{u')Mu)da ) dt 
Js \Jan N J 

For simplicity, let da m = (m.v)da. If we observe that E'(t) = — I g(u')u'da m < 0, we get, for a 

Jan N 

constant C > independent of E(0) if p = 1, 

2d f E^dt < CE(S) + f E 2 ^ ( f {u'f - \Vu\ 2 - g(u')Muda^\ dt. 
Js Js \Jan N J 

Using the definition of Mu and Young inequality, we get, for any e > 0, 

2d J T E P -^dt<CE{S) + J T E^ 1 ^ (V ) 2 + (iHlL + rf2( " fe ) 9(u') 2 + eu 2 ^ da m )dt. 

Now, using Poincare inequality, we can choose e > such that 

e I (m.v)u 2 da < ^ [ \Vu\ 2 dx < dE, 
Jan N 2 Jq 

so we conclude 

d f E E ^dt<CE{S)+C f E*r(( ((u) 2 + g{u') 2 ) da m )dt. 
j s j s J an?? 

We split dn 7v to bound the last term of the above estimate 

<9ft]v = {x e dQ N ; \u'{x)\ > 1}, dn 2 N = {i£ dCl N ; \u'{x)\ < 1}. 
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Using Q and (0]), we get 

r(| ({u') 2 + g{u') 2 ) dtr m )dt < C I U u'g{u')dcr m )dt < CE{S), 



T , f i-T 

p-1 

E— 

s x Jdn\, Js y Jdn? 



where C neither depend on E(0) if p = 1. 

On the other hand, using ([3]), Q; Jensen inequality and boundedncss of m. one successively obtains 



p+i 



((u') 2 +g{u') 2 ) da m < C / (u' 5 ( U ')) 2/(p+1) da m <C[ u'g(u')da m < C{-Et)vW. 

dn% Jdn% \JdQ% J 

Hence, using Young inequality again, we get for every e > 

/ E 1 ^ 1 ( f ((u') 2 + g(u') 2 ) dcr m ) dt < [ {sE 1 ^ -C{e)E')dt < e [ E^ dt + C(e)E(S), 
Js \Jdn 2 N J Js Js 

so finally, we get, for some C(e) and C independent of E(0) if p = 1 

d E E i 1 dt<C(e)E(S)+eC E^dt. 
Js Js 

Choosing now eC < ~, theorems results from Proposition [S] ■ 

4 Examples and numerical results 
4.1 Examples 

We here consider the case when Q is a plane convex polygonal domain. The normal unit vector pointing 
outward of fl is piecewise constant and the nature of boundary conditions involved by the multiplier 
method can be determined on each edge, independently of other edges. 

Along each edge, vector v is constant and the boundary conditions arc defined by the sign of 

m(x).v(x) = (Re{x - x )).v(x) = (x - x a ).R^ e (is(x)) . 

Hence we build is, R^g{v) and we we can determine the sign of above coefficient with respect to the 
position of xq. To this end, we construct two straight lines, orthogonal with respect to R-g{v) so that 
each of them contains one vertex of the considered edge. 

This determines a belt and if xq belongs to this belt, we obtained mixed boundary conditions along this 
edge, if xq does not belong to this belt, then we get Dirichlet or Neumann boundary conditions along 
whole the edge (see Figure 6). 



/ v / 



/ / 
/ / 

D ' N D / ___N_ 



/ / 
/ / 
/ / 

Figure 6: Boundary conditions along some edge depending on the position of xq. 

Performing this method for a square, fl = (0, l) 2 , we show in Figure 7 the different cases of boundary 
conditions depending on the position of x$. Three main cases are considered 

7T 

1. < 9 < — : above belts controlling opposite edges have a non-empty intersection, which is a belt 
of positive thickness, 
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this intersection is a straight line, 



7T 7T 

3. — < 9 < — : the intersection is empty. 
4 2 



The case when 9 is negative can be easily deduced by symmetry. 

In the three above cases, there are four angular sectors (shaded areas in Figure 7) such that if xq belongs 
to one of them, then geometrical condition © is satisfied. 



• J -J t 

□ 



C 



□ \ 



Figure 7: Shape of boundary data with respect to xq (from left to right, cases 1,2,3). 
4.2 Numerical results 

We perform numerical experiments by considering the following case 

ft = (o,i) 2 , dn D = ({o}x [o,i])u([o,i] x{o}), m N = asi\an D , 

and using above vector field 

m(x) = Rg(x — xq) . 
We only consider the case of a linear feedback. Let us write down the problem. 

' u" - Au = , infixi;, 

u = , on dU D x R* + , 

d v u = — (m.v) u' , on 8£In X , 

u(0) = Wo , in fi , 

u (0) = u\ , in . 

We will investigate cases when 9 varies in [0, arctan(2)]. A particular case is given in Figure 8. 





a 


n 













on. 



Figure 8: When xq belongs to Vg 1 geometrical condition © is satisfied at a. 

Our aim here is to study numerically the variations of the speed of stabilization with respect to the 
position of Xq and the value of 9. 



1G 
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To this end, we have built a finite differences scheme (in space). This leads to a linear second order 
differential equation 

U" + BU' + KU = , (18) 

where B is the feedback matrix and —K is the discretized Laplace operator. 
This differential equation can be rewritten as follows 




and the energy function is approximated by -((U,KU) + \\U'\\ 2 ) = -{\\V\\ 2 + \\U'\\ 2 ) with V = K 1/2 U. 

The decreasing rate is given by the highest eigenvalue of above matrix. Results of our computations are 
shown in Figure 9 where we built the decreasing rate as a function depending on 9 and the position of 
Xq represented by the abscissa A along T>g. 



0.07 






theta lambda 



Figure 9: Dependancc of the decreasing rate with respect to 6, A. 

It can be observed that in this case, the decreasing rate is increasing with 9 and the best position for xq 
is the origin of half-line Vg . 

References 

[1] Brezis, H., 1973, Operateurs maximaux monotones et semi-groupes de contractions dans les espaces 
de Hilbert, Math, studies 5, North Holland. 

[2] Banasiak, J., Roach, G.-F., 1989, On mixed boundary value problems of Dirichlet oblique- 
derivative type in plane domains with piecewise differentiable boundary. J. Diff. Equations, 79, no 
1, 111-131. 

[3] Bardos, C, Lebeau, G., Rauch, J., 1992, Sharp sufficient conditions for the observation, control 
and stabilization of waves from the boundary. SIAM J. Control Optim., 30, no 5, 1024-1065. 



Cornilleau, Loheac, Osses 



17 



[4] Bey, R., Loheac, J. -P., MOUSSAOUI, M., 1999, Singularities of the solution of a mixed problem 
for a general second order elliptic equation and boundary stabilization of the wave equation. J. Math, 
pures et appli., 78, 1043-1067. 

[5] Conrad, F., Rao, B., 1993, Decay of solutions of the wave equation in a star-shaped domain with 
non linear boundary feedback. Asymptotic Analysis, 7, no 1, 159-177. 

[6] Grisvard, P., 1985, Elliptic problems in nonsmooth domains. Pitman, London. 

[7] Grisvard, P., 1989, Controlabilite exacte des solutions de l'equation des ondes en presence de 
singularites. J. Math, pures et appli., 68, 215-259. 

[8] Kozlov, V. A., Maz'ya, V. G., Rossmann, J., 1997, Elliptic Boundary Value Problems in Do- 
mains with Point Singularities, AMS, Providence. 

[9] Komornik, V., 1994, Exact controllability and stabilization; the multiplier method. Masson-John 
Wiley, Paris. 

[10] Komornik, V., Zuazua, E., 1990, A direct method for the boundary stabilization of the wave 
equation. J. Math, pures et appi, 69, 33-54. 

[11] Lions, J.-L., 1988, Controlabilite exacte, stabilisation et perturbations de systcmcs distribucs. 
RMA 8, Masson, Paris. 

[12] Ho, L.F., 1986, Obscrvabilite fronticre de l'equation des ondes. C. R. Acad. Sci. Paris, Ser. I Math. 
302, 443-446. 

[13] MOUSSAOUI, M., 1996, Singularites des solutions du probleme mele, controlabilite exacte et stabilisa- 
tion frontiere. ES AIM Proceedings, Elasticite, Vis co elasticite et Contrdle optimal, Huitiemes Entretiens 
du Centre Jacques Cartier, 157-168. 

[14] Necas, J., 1967, Les methodes directes en theorie des equations clliptiques. Masson, Paris. 

[15] OSSES, A., 1998, Une nouvelle famille de multiplicateurs et applications a la controlabilite exacte 
de l'equation des ondes. C. R. Acad. Sci. Paris, 326, Serie I, 1099-1104. 

[16] OSSES, A., 2001, A rotated multiplier applied to the controllability of waves, elasticity and tangential 
Stokes control. SI AM J. Control Optim., 40, no 3, 777-800. 

[17] Rellich F. , 1940, Darstellung der Eigenwerte von Au+Xu durch ein Randintegral. Math. Zeitschrift, 
46, 635-636. 

[18] Shamir, E., 1968, Regularity of mixed second order elliptic problems. Israel Math. Journal, 6, 
150-168. 



